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NOTATION 

Complex  wave  amplitude 
Depth  of  submergence  of  wedge 
=  dN/U 

Function  defined  by  Equation  [49] 

Function  defined  by  Equation  [6o] 

Group  velocity 
Heaviside  function 

Components  of  group  velocity  In  Cartesian 
coordinates 

Components  of  group  velocity  In  plane  polar 
coordinates 

Components  of  group  velocity  In  cylindrical 
polar  coordinates 

Radial  component  of  wave  number  in  plane  polar 
coordinates 

Radial  component  of  wave  -"umber  in  cylindrical 
polar  coordinates 

Half-length  of  wedge 

Source  strength 

Vaisala  frequency 

Zero  or  positive  integer 

Source  distributions 

Fourier  transform  of  q,  qi  and  qa  respectively 
Radial  coordinate  in  cylindrical  polar  coordinates 
Radial  coordinate  in  plane  polar  coordinates 
Parameter  in  Laplace  transform 
Time 
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U,  ui,  ua 
Xj  y,  z 

y 

y* 


a,  /3,  7 


9 

e 

T 

0) 


Horizontal  velocity  components 
Cartesian  coordinates 

Also  an  axial  coordinate  In  cylindrical  polar 
coordinates 

=  yN/U 

Cartesian  components  of  wave  number 

Polar  angle  of  wave  number  plane  In  plane  polar 
coordinates 

Polar  angle  of  wave  number  space  In  cylindrical 
polar  coordinates 

Half  of  wedge  angle 

Polar  angle  In  cylindrical  polar  coordinates 
Polar  angle  In  plane  polar  coordinates 
Frequency 
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INTRODUCTION 

In  a  recent  paper,  Wu  and  Mel  (196?)  discussed  the  problem 
of  gravity  waves  generated  by  a  submerged  two-dimensional  dis¬ 
turbance  moving  horizontally  in  a  stratified  fluid  with  a  free 
surface  using  a  linearized  theory.  In  addition  to  the  usual 
free  surface  wave  mode.  Internal  waves  that  behave  asymptotically 
like  outgoing  cylindrical  waves  were  found  on  the  downstream  side 
of  the  disturbance.  Experimental  evidences  (Long  1955j  Ylh 
1959)>  however,  indicate  that  unattenuated  waves  (blocking)  can 
also  exist  upstream  as  well  as  downstream  of  the  disturbance. 
Since  Wu  and  Mel  have  only  examined  the  vertical  component  of  the 
perturbation  velocity  their  analyses  do  not  rule  out  the  possi¬ 
bility  that  a  linear  theory  can  account  for  the  existence  of  the 
blocking  phenomenon;  the  reason  for  this  is  that  the  blocking 
perturbations  would  consist  essentially  of  horizontal  motions. 

This  report  attempts  to  resolve  this  issue  within  the  frame¬ 
work  of  a  linearized  theory  incorporating  the  Bousslnesq  and 
Oseen  approximations.  The  fluid  will  be  taken  to  be  unbounded 
and  possess  a  constant  Vaisala  frequency.  The  concept  of  group 
velocity  (Lighthill  i960,  1967)  will  be  used  to  clarify  the 
physical  basis  of  the  solution.  In  addition  to  the  two-dimen¬ 
sional  problem,  three-dimensional  effects  on  blocking  will  be 
considered.  In  the  last  section  the  internal  waves  generated 
by  a  two-dimensional  wedge  moving  horizontally  in  a  stratified 
fluid  beneath  a  surface  of  large  density  discontinuity  is  con¬ 
sidered. 
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PORMULATION  OP  THE  PROBLEM 

Consider  the  problem  of  a  submerged  disturbance  moving  with 
constant  horizontal  velocity  of  magnitude  U  in  an  infinite  ex¬ 
panse  of  density-stratified,  incompressible  and  invlscid  fluid 
under  gravity.  The  Valsala  frequency  N  of  the  fluid  is  assumed 
to  be  constant.  Choose  a  Cartesian  coordinate  system  (x,  y)  with 
origin  at  the  center  of  the  disturbance  and  oriented  in  such  a 
manner  that  gravity  points  in  the  negative  z  direction  and  the 
fluid  appears  to  be  moving  in  the  positive  direction. 

Let  us  assume  that  the  disturbance  can  be  represented  by  a 
source  distribution  q(x,y).  A  large  class  of  moving  solid  bodies 
can  thus  be  represented,  although  the  correspondence  is  not 
simple  since  the  body  shape  depends  in  a  complex  way  upon  the 
stratification  of  the  ambient  fluid  and  velocity. 

Within  the  Boussinesq  and  Oseen  approximations  the  linearized 
equation  governing  the  horizontal  component  of  the  perturbation 
velocity,  u,  is 


{11 

where  t  is  the  time.  We  focus  our  attention  on  this  component 
of  the  velocity  because  we  expect  the  perturbations  far  upstream 
and  downstream  of  the  disturbance,  if  there  are  any,  to  be  es¬ 
sentially  horizontal.  For  a  three-dimensional  disturbance  the 
corresponding  equation  would  be 
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A 

ht 


+ 


1  ^  \  n 

,2 

u  = 

dx®  ^z^  J 

[2] 


where  z  is  the  third  space  coordinate. 


It  is  well  known  that  problems  of  the  sort  being  formulated 
are  not  determinate  and  it  is  necessary  to  have  a  way  to  select 
that  particular  solution  corresponding  to  physical  reality.  Fol¬ 
lowing  Wu  and  Mel  (1957)  we  consider  the  corresponding  transient 
problem  in  which  the  forcing  function  is  switched  on  at  t  «  0 
and  take  the  large  time  limit  of  the  resulting  time-dependent 
solution.  The  switching  operation  is  accomplished  by  multip''ying 
the  source  term  by  the  Heaviside  function  H(t). 


FEATURES  OF  THE  WAVE  PATTERN 

Without  Integrating  (ll  directly  we  can  obtain  certain  in¬ 
formation  about  its  associated  regular,  steady  wave  pattern 
using  the  concept  of  group  velocity.  The  regular  wave  pattern 
may  be  considered  to  be  composed  of  plane  wave  components  and, 
because  the  pattern  is  steady,  only  those  plane  waves  that  are 
stationary  with  respect  to  the  disturbance  qualify  to  be  in¬ 
cluded.  Packets  of  these  waves  are  continuously  generated 
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through  the  Interaction  of  the  disturbance  with  the  free  stream. 
These  wave  packets  are  then  swept  downstream  and,  at  the  same 
time,  they  propagate  with  their  own  group  velocities. 

Let  us,  therefore,  consider  the  behavior  of  plane  waves  in 
a  stratified  fluid  stream  to  select  those  which  qualify  to  be  in¬ 
cluded  in  the  regular  wave  pattern.  Two  dimensional  plane  waves 
have  the  form 


i(ax  +  |3y  -  (ut) 

u  =  A  e  (31 

where  the  components  of  the  wave  number  vector  (a,/5)  and  the 
frequency  oi  are  real  numbers,  and  the  amplitude  A  Is  a  complex 
number.  They  must  satisfy  the  associated  homogeneous  Equation  of 
[ll.  This  leads  to  the  following  dispersion  relation: 

(tu  -  Ua)*(a*  +  i3*)  -  -  0  (4) 


which  has  two  branches,  viz. 


u)x  »  Uo - j:"  (5) 


Nq 


u>a  ■  Ua  + 


(6) 
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All  wave  numbers  corresponding  to  osi  =0  and  cua  =  0  represent 
stationary  plane  waves,  they  are 


a  =  0  for  0)1  =  0  and  0)2  =  0  [7] 


and  are  shown  in  Figure  1  where  a  double  line  is  used  to  Indicate 
the  double  root  at  a  =  0.  ""he  regular  wave  pattern  will  be  com¬ 
posed  solely  of  these  stationary  plane  waves. 

To  obtain  further  information  about  the  wave  pattern  we 
use  the  concept  of  group  velocity.  Let  us  consider  first  the 
system  of  plane  waves  corresponding  to  the  double  root  at  a  =  0. 
They  have  wave  numbers  pointing  in  the  vertical  direction  and 
their  group  velocities  are  given  by 

,  0  =  0  for  0)1  *  0  [9) 

0  »  +  —  0  .=  0  for  0)2  =  0  [10] 

Where  and  denote  components  in  the  x  and  y  directions 
respectively.  I'ackets  of  these  waves,  therefore,  i^ropagate  un¬ 
attenuated  in  the  horizontal  direction.  Since  the  wave  packets 
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are  swept  downstream  with  velocity  (U, O),  waves  corresponding  to 
the  branch  cua  =  0  always  trail  behind  the  disturbance  while  those 
corresponding  to  cui  =0  appear  ahead  of  the  disturbance  if 

jjSj  <  ^  and  behind  if  j |  ^  .  When  [^j  =  ^  the  group 

velocity  is  equal  in  magnitude  but  opposite  in  direction  to  the 
free  stream  velocity.  These  wave  packets,  therefore,  do  not  move 
away  from  the  disturbance  after  their  generation  and  thus  their 
amplitudes  grow  with  time.  The  linear  theory  is  not  adequate  for 
describing  the  behavior  of  these  wave  components. 

Next  we  consider  the  system  of  plane  waves  corresponding  to 
[6]  —  the  cylindrical  system.  It  is  convenient  to  introduce 
polar  coordinates  defined  as 


a  =  k  cos  0, 
X  =  r  cos  <t>. 


J3  =  k  sin  9 
y  =  r  sin  <)> 


[11] 


In  polar  coordinates  the  group  velocities  (G  ,G.  )  of  the  com- 

1*  9 

ponents  of  this  system  of  waves  are  given  by 


G  =0,  G«  =  U  sin  e  [12] 

r  0 

Figure  2  shows  the  vector  sum  of  the  group  velocity  and  the  free 
stream  velocity  for  various  components  of  this  system.  It  is 
seen  that  they  all  point  in  the  direction  downstream  of  the  dis¬ 
turbance.  Thus,  this  system  of  waves  always  trails  behind  the 
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dlsturbance.  When  9  =  ±  ~  the  group  velocity  is  exactly  equal 
In  magnitude  but  opposite  in  direction  to  the  free  stream  ve¬ 
locity  so  that  again  the  linear  theory  is  invalidated. 


The  resulting  wave,  of  course,  also  depends  on  the  excitation 
produced  by  the  disturbance.  For  the  linear  theory  to  be  valid 
we  need  to  restrict  ourselves  to  those  source  distributions  whose 
Fourier  transforms  vanishes  at  the  critical  points 


a  =  0, 


Let 


[13] 


q  =  qi  +  q2 


and  using  a  bar  over  a  function  to  denote  its  I^ourler  transform 


^  we  have 

q  =  qi  +  Qs 

[15 

The  above  requirement  is  satisfied  if  qi  vanishes  at  a  =  0  and 

5  a 

qa  at  a  .+  .S  -  -r  =  0  because  the  sum  q  vanishes  at  the 

U 


critical  points;  its  values  at. any  other  point,  however,  is 
arbitrary. 


This  limitation  of  the  theory  is  entirely  due  to  finite 
amplitude  effects  and  has  nothing  to  do  with  the  Bousslnesq  ap¬ 
proximation.  Following  an  argument  similar  to  the  above,  one  can 
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easlly  show  that  elimination  of  the  Boussinesq  approximation 
fails  to  remove  this  limitation. 

THE  REGULAR  WAVE  PATTERN 

In  the  previous  section  we  have  obtained  a  broad  picture  of 
the  regular  wave  pattern,  here  we  will  derive  it  analytically. 
Referring  to  [1]  let  us  split  u  into  two  components  ui  and  U2 
such  that 


and 


qiH{t) 


[16] 


q2H(t) 


[17] 
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It  is  noted  that  we  have  split  the  total  source  distribution 
into  two  components  qi  and  qz  as  discussed  in  the  previous  sec¬ 
tion  where  qi,  the  Fourier  transform  of  qi  vanishes  at  a  =  0  and 

q2  at  a®  +  jS®  -  —  =  0. 

Using  the  method  of  Pourier-Laplace  transforms  we  obtain 
the  following  integral  representations  of  the  steady  state  solu¬ 
tions  of  [l6]  and  [1?]: 


U2 


a[ (s+iaU)^+  N^]  qi 
(s+iaU)^(a®+  i3®)+ 


i(ax+  i3y) 
e 


[18] 


ar(.-HaU)^  f  N»1  a% 

(s+iaU)®(a^+  i3^)  + 


[191 


Let  us  first  consider  ui  and  rewrite  the  integral  in  polar 
coordinate  form  as  follows: 


7r/2  « 

^m  fdefdk  cos  e)^-t  N^]qi(k  cos  6.k  sin  e).. 

2w^  s-»  0  J  J  (s  +  ikU  cos  9)®  +  N®cos®  6 

0  0 


ikr  cos  (0-^)  ikr  cos  (6+<>) 
e  +  e 


[20] 
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using  the  method  of  stationary  phase  we  obtain  the  following 
asymptotic  form  for  ui ; 
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From  [26]  It  Is  noted  that  disturbances  whose  Fourier 
transforms  vanish  at  a  =  0  cannot  excite  the  fore  and  aft  unat¬ 
tenuated  wave  system  and^  therefore,  no  blocking  can  occur.  These 
source  distributions  satisfy  the  condition 

00 

/  q(x,y)  dx  =  0  [27] 


that  is,  for  all  y  the  sources  are  exactly  balanced  by  the  sinks. 
In  particular,  a  source  distribution  that  is  antisymmetric  in  x 
would  satisfy  this  condition. 

It  is  also  noted  that  the  unattenuated  waves  ahead  of  the 
disturbance  have  small  wave  numbers.  Thus  a  disturbance  whose 
vertical  extent  i  is  small  compared  with  U/^I  can  excite  only 
very  weak  forward  u.. attenuated  waves  because  the  Fourier  trans¬ 
form  of  the  disturbance  will  be  small  in  the  low  wave  number 
range. 

By  integrating  ua  over  vertical  sections  upstream  and  down¬ 
stream  of  the  disturbance  it  is  noted  that  one  half  of  the 
total  discharge  from  the  source  distribution  flows  upstream  and 
the  other  half  flows  downstream. 

THREE-DIMENSIONAL  EFFECTS 

We  have  seen  above  that  fore  and  aft  unattenuated  waves  can 
be  excited  by  two-dimensional  disturbance?  moving  horizontally 
in  a  stratified  fluid.  For  a  disturbance  having  a  finite  span- 
wise  extension,  however,  this  wave  system  attenuates  due  to 
lateral  spreading  and  therefore,  no  blocking  can  occur. 


HYDRONAUTICS,  Incorporated 


The  equation  corresponding  to  [4]  for  the  three-dimensional 


case  is 


(a3-Ua)^(a^+  j3^+  y^)  -  N®  (a^+  7®)  =  0 


which  also  has  two  branches,  viz. 


Oil  =  Ua  -  N 


a^+  7= 


a^+  /3®+  7® 


[29] 


0)2  =  Ua  +  N 


a2+  7® 


a®+  ^®+  7® 


[30] 


In  polar  cylindrical  coordinates  defined  as 


a  «  K  cos  0;  7  “  K  sin  0;  /5  »  ^ 


X  =  R  cos  z  ■  R  sin  <>;  y  »  y 


[31] 


they  are 


0)1  ■  UK  cos  0  -  N 


K®  +  ^* 


[32] 


^K*  + 


t33l 


0)2  ■  UK  cos  0  +  N 
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Therefore,  the  stationary  plane  waves  are 

K  =  0  for  (Oi  =  0  and  0)2  =  0  [3^1 

~^K^+  cos  6  =  ^  for  (Oi  =  0  cos  0  >  0  [35] 

cos  0  =  -  ^  for  u)2  =  0  cos  0  <  0  [36] 

The  system  of  plane  waves  given  by  [3*^1  corresponds  to  the  wave 
system  that  causes  blocking  in  the  two-dimensional  case.  Our 
primary  interest  is  to  see  the  three-dimensional  effects  on  the 
blocking  phenomenon  and  thus  we  shall  concentrate  our  attention 
on  Phis  system. 

Members  of  this  system  have  z  »  constant  planes  as  surfaces 
of  constant  phase  and  the  components  In  cylindrical  coordinates 
of  their  group  velocities  are 

Op  ^  ,  0  «  0,  0  ■  0  for  0)1-0  [37) 

R  p  0  y 

Op  -  +  —  ,  0-  -  0,  0„  -  0  for  0)8  -  0  [38I 

R  ^  e  y 

It  is  seen  that  those  corresponding  to  a)j  -  0  have  energy  flow¬ 
ing  radially  Inwards  and  are,  therefore,  not  excited  by  a  lo¬ 
calized  disturbance. 
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The  wave  packets,  while  propagating  with  their  own  group 
velocities  are,  at  the  same  time,  being  washed  downstream.  Fig¬ 
ure  3  shows  the  plane  view  of  the  position  of  wave  packets  of 
different  \^\  one  unit  of  time  after  *.aeir  generation  by  the  dis¬ 
turbance  which  Is  localized  at  the  origin  of  the  coordinate  sys¬ 
tem.  Wave  packets  with  |j3|  >  N/U  trail  behind  the  disturbance 

and  are  restricted  to  a  sector  with  a  half  Included  angle  of 

- 1  N  I  I  N 

sin  while  packets  with  |/3|  appears  in  front  as  well 

as  behind  the  body. 

THE  WAVES  GENERATED  BY  A  TWO-DIMHJSIONAL  WEDGE 

In  this  last  section  we  would  like  to  apply  some  of  the 
theoretical  results  developed  In  the  previous  sections  to  a 
specific  problem. 

As  we  have  seen,  the  Fourier  transform  of  the  source  dis¬ 
tribution  that  represents  the  dlstur’^an-'e  must  vanish  at  the 
critical  points  given  by  Equation  C 1.  ]  for,  otherwise,  waves  that 
have  group  velocltl'  equal  but  qpposlte  to  the  main  stream 
would  be  excited.  The  linear  theory  Is  not  adequate  for  treating 
these  particular  wave  components.  One  of  the  source  distribu¬ 
tions  that  satisfies  this  requirement  consists  of  two  line  sources 
of  constant  strength  m  arranged  as  ,thown  In  Figure  4  where  21 
denote  the  length  of  the  source  lines  and  2d  Is  the  distance 
separating  them.  It  will  later  be  shown  that  d  has  to  be  re¬ 
lated  to  U  and  N  by 
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d  =  7r(n  +  U/N  l39] 

where  n  is  zero  or  any  positive  integer.  Taking 

m  -  213  tan  t  [40] 

where  t  Is  supposed  to  be  small,  then  according  to  the  thin  body 
theory,  the  body  corresponding  to  each  of  the  source  lines  is  a 
wedge  of  angle  2t  and  length  21  folic  1  by  a  stern  of  thickness 
2i  tan  T.  If  both  t  and  I  tan  t  are  small  then, aside  from  the 
vicinity  of  the  stagnation  point  and  the  turning  point  at  the 
shoulders, the  perturbations  are  small  and  we  should  expect  the 
linear  theory  to  be  reasonably  adequate.  Furthermore,  from  the 
governing  equations  and  the  symmetry  of  the  disturbance,! t  is 
clear  that  the  flow  is  symmetric  about  the  x-axis  and  the  flow  in 
the  lower  half  plane  may  be  considered  to  be  produced  by  a  wedge 
moving  horizontally  beneath  a  wall  or  a  surface  of  large  density 
discontinuity. 

The  Fourier  transform  of  this  source  distribution  Is 

q  »  4m  cos  3d  [41] 


At  a  ■  0 


q(0,3)  •  4ml  cos  3d 


(42] 
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therefore  d  has  to  satisfy  Equation  [39]  in  order  that  q{0,^) 
vanishes  at  the  critical  points. 

q  can  now  be  written  as  the  sum  of  qi  and  qa  so  that  qi 
vanishes  at  a  =  0  and  qa  at  -  N^/U®  =  0  as  follows 


qi 


AmU 

N  sin  9 


cos 


m 

U 


cos  0 


sin 


.  U 


sin 


6 


[45] 


qa  =  4m  cos  /3d 


sin  all 


4raU 


a 


N  sin  6 


cos 


Nd  . 
—  cos  9 


sJn 


'M. 

,  U 


sin  0 


[46] 

The  outgoing  wave  system  trailing  behind  the  wedge  is  ob¬ 
tained  by  substituting  equation  [45]  into  equation  [24]  to  give 


Ui  ~  H(x)2m 


tan  <t> 


f  Nd  ^  1 

[Nr  TT 

cos  1 

j-  cos 

j  sin  j 

TT-sln  <t)  1 

i 

cos 

|u  -  *) 

[4?] 

We  shall  not  get  into  the  details  of  this  wave  system  beyond 
mentioning  the  fact  that  the  wave  is  not  singular  at  <l>  *  ±  7r/2, 
l.e.,  along  the  y  axis,  in  spite  of  the  factor  tan  because 
cos  (Nd  cos  <I)/U)  vanishes  there. 
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where 

F2(y*,d*)  =  -  sin  d*  cos  y*[sgn(d*  +  y*)  +  sgn  (d*  -  y*)] 

[51] 

In  Figure  6,  the  expression  (Pi  -  irFs)  Is  plotted  against  y* 
for  d*  given  by  Equation  [50]. 

It  is  noted  that  for  the  case  where  the  fore  and  aft  unat- 
'tenuated  perturbation  profiles  depend  only  on  the  volume  of  fluid 
displaced  by  the  wedge  (l.e.  ml)  but  not  on  the  shape  of  the 
wedge,  we  can  arbitrarily  vary  m  and  I  respectively  so  long  as  we 
keep  ml  constant.  The  cylindrical  wave  system,  however,  would 
not  stay  the  same.  - 
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FIGURE  1  •  WAVE  NUMBERS  OF  STATIONARY  PUNE  WAVES 
IN  A  UNIFORM  STREAM 
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FIGURE  2  -  VELOCITY  OF  PROPAGATION  OF  WAVE  PACKETS  OF 
CYLINDRICAL  SYSTEM 
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FIGURE  3  -  POSITION  OF  WAVE  PACKETS  OF  DIFFERENT  181  ONE  UNIT  OF 
TIME  AFTER  GENERATION 
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FIGURE  4  -  SOURCE  DISTRIBUTION  REPRESENTING  WEDGE  BENEATH  SURFACE  OF 
LARGE  DENSITY  DISCONTINUITY 


HGURE  5  -  VaOQTY  PERTURBATION  UPSTREAM  OF  WEDGE 
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nCURE  6  •  VELOQTY  PERTURBATION  DOWNSTREAM  OF  WEDGE 
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